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METHOD AND DEVICE FOR THE ESTIMATION OF A PROPAGATION 
CHANNEL FROM ITS STATISTICS 

BACKGROUND OF THE INVENTION 
1. Field of the Invention 
5 The present invention relates to a method for the estimation of a 

propagation channel from its statistics. 

It can be applied in the third-generation mobile telephony networks 
commonly known as UMTS networks. 

It can also be applied to multiple-sensor receivers used for 
10 functions of metrology, capacity augmentation, or again for radio-monitoring 
, A functions. 

□ 2. Description of the Prior Art 

)% In a transmission system, especially one using radio waves, a 

y transmitter sends out a signal in a transmission channel to a receiver. The 

J* 15 signal that is sent undergoes amplitude and phase fluctuations in the 
tl transmission channel. The signal received by the receiver consists of time- 

shifted and modified copies of the signal sent. The fluctuations of the signal 
and the shifts generate interference known by those skilled in the art as 
ry intersymbol interference. This interference arises especially from the law of 

H 20 modulation used for the transmission and also from multipath propagation in 
fjj the channel. 

The received signal generally results from a large number of 
reflections in the channel, the different paths taken by the transmitted signal 
leading to various delays in the receivers. Thus, the impulse response of the 
25 channel represents all the fluctuations undergone by the transmitted signal. 

The forthcoming arrival of UMTS networks is obliging equipment 
manufacturers and suppliers to adapt metrology, capacity-augmentation and 
radio-monitoring tool3 to this new standard. The performance characteristics 
of these tools rely partly on the estimation of the effect of the propagation 
30 medium (the radio channel) on the signals sent by the different entities of the 
network (base station, mobile units etc). This operation is known as 
propagation channel estimation. It can be used especially to counter the 
effect of the propagation to improve the quality of the signal received (by 
equalization) or to bring out information on the propagation medium (namely 
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the directions of arrival and the path delays) in order to implement spatial 
processing operations such as direction-finding. 

Modelling 

Signal sent 

i The composition of the UMTS signal sent by the base station to 

the mobile units (in the downlink) may be modelled for example according to 
the scheme shown in figure 1. 

The signal s(t) consists of several frames, each frame having a 
fixed duration and comprising a given number of slots. For example a10ms 
frame comprises 15 slots, 

Figure 2 gives a diagrammatic view of a modulator for the 
downlink (from the base station to mobile units). 

The Q-PSK (Quadrature Phase Shift keying ) symbols, referenced 

bq(0) bq(Ns-1), intended for a user q, are first of all multiplied by a power 

factor jjq, with Ns being the number of symbols sent to this user. Each 
symbol is then modulated by a sequence known as a spreading sequence, 
referenced cq with a value of ±1 , and a size Nq (spreading factor). Thus, the 
symbol b q (l) is used to form the sequence // q b q (l) c q (0), jjq b q (l) Cq(N q - 
1 ), where I is the index of the symbol. 

The sequences cq are orthogonal so that: 
2c,(«)cv0.)=^ if p = q 

Else it is equal to 0 

The sequences thus formed are then multiplied term by term by a 
sequence of symbols , (±1+i), called a scrambling code s. This sequence is 
periodic, with a frame period, and is built so as to simulate a random signal. 

The resulting signal to be sent is therefore written as follows: 
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with s(n) being the scrambling at a given point in time n 
where LJ designates the "integer part" operator and [ ] designates 
the "modulo" operator. This can also be written as: (1) 



with ; 

L"«J (2) 
This signal goes to a transmission filter and is then broadcast by 
5 the base station to the mobile units. 

Signal received bv a mobile unit 

The signal is received on a network of sensors after crossing the 
radio channel. After sampling, it can take the form: 

x(«) = £;K«-*)hC*)+b(«) 
t *=» (3) 

53 io * where the boid characters designate vectors, and 

|S * L corresponds to the spread of the channel expressed in number 

W of chips, the indices n and k correspond to chips . 

jfy * x(n) represents: 

S * when the sampling is done at the chip rate , x(n) is the vector 

q is of the signals received at the instant n on each sensor, 

FU * when over-sampling is done at a rate corresponding to one 

:H chip/2, x(n) is the vector of the signals received at the instants 

p n and n+chip/2. 

ru * h(k) is the multi-sensor channel, and 

20 * b(n) is an additive noise combining the interference phenomena 

coming from the other base stations and the thermal noise. 

Signal received by a station 

Figure 3 gives a diagrammatic view of one possibility of modelling 
the uplink (from the mobile units to base station ) described in detail here 
25 below in the context of the method according to the invention. 
Known methods 

There are many known methods for estimating the response of a 
propagation channel. 

For example, one classic method proceeds by correlation of the 
30 received signal with shifted versions of a known learning sequence of the 
receiver. Of this sequence, orrty Ihe steps needed for the understanding of 
the invention are recalled here below. 
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For the user 0, for example, the channel impulse response is 
estimated, in a classic manner, by correlating several shifted versions of the 
received signal x(n+k) (one version being shifted by the sample k) by the 
learning sequence; 



(4) 

where P represents the number of learning symbols, and 
NoP the number of learning chips, or the number of pilot chips 
By combining the equations (1), (2), (3) and (4), the detailed 
expression of the estimated channel corresponds to (5) : 

15 

The properties of the sequences dq(n) are such that (6) : 
- . L , sinon 

20 Tl^P 



Thus, the estimate of the channel impulse response is expressed 
as a function especially of the response of the channel h (k) and of the three 
25 terms representing existing interference (7) : 

1 tf,p~: 



* where Q is the number of users of the propagation channel, 

* I, k are indices corresponding to chips, and 



• n is the index of the instant considered 

Three types of interference disturbing the estimation are thus 

identified: 

• the first term corresponds to the self-correlation for the 
5 sequence d Q , 

• the second term corresponds to the interference from the other 
users, and 

• the third term corresponds to the contribution of the external 
noise and the thermal noise. 

10 This technique is efficient when the (earning sequence is long and 

when the propagation channel does not change or undergoes little change in 
K time. In the case of a fast variation of the propagation channel, it becomes 

O necessary to estimate it on fairly short periods of time. 

fi Another technique, known as the "least error squares technique" 

£ 15 improves the above method by removing the need for self-correlations of the 
learning sequence. 

The patent FR 2 762 164 discloses a method to estimate the 
O impulse response of the transmission channel. This method uses the 

jj estimation of the space-time covariance matrix r of the impulse response of 

S! 20 the channel. The method considers that the estimation noise is white, with a 
fjj power B, and independent of the channel. In this case, the matrix r is 

estimated by f* = A-Bl with A being the matrix of space-time covariance of 
the estimated channel. If B is unknown, it can be estimated by the smallest 
eigenvalue of the estimated matrix A. It can also be fixed at a threshold 
25 value. Such an assumption is perfectly suited to a system in which a 
propagation channel is dedicated to a user, for example the GSM system. 
Furthermore, the proposed technique relates to single-sensor receivers. 

However, a method of this kind is no longer suited to signals 
including several users, such as the UMTS signals where the noise is neither 
30 white nor independent of the channel, especially owing to multiple-user 
interference. 

An object of the present invention is a method to estimate the 
propagation channel from its statistics, which are themselves estimated by 
expressing especially the noise from the matrix of empirical correlation of the 
35 observations. 



SUMMARY OF THE INVENTION 

The invention relates to a method to estimate the impulse 
response h of a propagation channel in a system comprising at least one or 
more sensors. The method comprises at least one step for estimating the 
5 statistics of the additive noise resulting from the interference and from the 
thermal noise on the basis of the statistics of the received signal. 

The step for estimating the covariance matrix of the noise is 
carried out, for example, on the basis of the empirical covariance matrix of 
the observations Rx and the number of pilot chips of a learning sequence 
10 transmitted with the signal, the noise matrix being expressed in the form 

J4 The method is used for example to estimate the impulse response 

Ul of a propagation channel in the UMTS field for uplinks and/or downlinks 

p between a base station and one or more mobile units, 

%i 15 The invention also relates to a reception device adapted to 

% s estimating the impulse response of a propagation channel. The device 

ry comprises one or more sensors for the reception of the signal, a means to 

rU sample the received signal, a means adapted to estimating the noise from 

the statistics of the channel. 
fU 20 The receiver according to the invention may comprise a means 

adapted to estimating the noise from the empirical covariance matrix of the 
observations R * and from the number of pilot chips of a learning signal 
transmitted with the signal, the matrix of the noise being expressed in the 
following form: 
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According to one alternative embodiment, it comprises for 
example a means adapted to determining the impulse response of the 
channel in the form: 

30 

H ra 4 i~ TH w Hf- — RJA-'H, 
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The receiver is used, for example, in the field of UMTS. 

The invention has ihe advantage especially of estimating a 
propagation channel comprising several users, for which the noise is neither 
white nor independent. 
5 BRIEF DESCRIPTION OF THE DRAWINGS 

Other features and advantages of the invention shall appear more 
clearly from the following description in which exemplary modes of 
implementation are proposed by way of an illustration, with reference to the 
appended figures, of which: 
10 * Figure 1 shows the composition of the UMTS signal, 

• Figure 2 shows the modelling of the signal for a downlink, 

- Figure 3 is a schematic drawing of a first alternative 
implementation of the invention, 

• Figure 4 shows an exemplary signal for the uplink (between a 
15 mobile station and a base station), 

• Figure 5 shows a second alternative mode of implementation 
of the invention, 

• Figure 6 shows comparative results obtained by using different 
methods of estimation 

20 In order to provide for a clear understanding of the object of the 

present invention, the following description, given by way of an illustration 
that in no way restricts the scope of the invention, is applied in the UMTS 
field to estimate a propagation channel especially on the basis of its 
statistics. 

25 The term " statistics" designates, for example, the space-time 

covariance of the impulse response of the channel. 

This system makes use of signals such as those described 
respectively with reference to figure 1 for a downlink, namely the signal sent 
from the base station to the mobile units, and with reference to figure 5 for an 

30 uplink corresponding to the signal sent from a mobile unit to the base station. 

By applying a realistic statistical model to the propagation channel, 
namely a model based on the assumption that the propagation channel can 
be likened to a finite impulse response filter, it is possible to improve the 
existing techniques of channel estimation, in taking account of the statistics 

35 of the channel. 
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Modelling of the redio channel 

The radio channel can be modelled as a sum of paths, each path 
being characterized by an arrival time, an attenuation and a phase shift. The 
modelling considered in the context of the invention uses the assumption 
5 according to which the attenuation and the phase shift of a path change 
rapidly as compared with the changes in the arrival time of the different 
paths. The channel may then be represented as a stationary random 
process. Thus, it can be said that each implementation of the channel is 
governed by the same statistics (the term "implementation" herein designates 
10 the value taken by the channel on a given slot, i.e. H m ). 

Thus, for a given slot with an index m , where m varies from 1 to 
M, for example, the vector of the coefficients of the sampled channels 
corresponding to each sensor of the receiver is referenced (8); 

s Hm=[h m (0) H h m {1) H ... h m (L) H J H 

Where the exponent H corresponds to the Hermitian transpose, L 
corresponds to the number of samples, and m is the index of the slot of the 
signal received. 

It is assumed that H m remains constant throughout the duration of 
the slot, but that for nwn', Hm^Hm' 

Similarly, the channel is estimated at each slot m, and its estimate 
is referenced The invention uses for example the method of successive 
correlations known to those skilled in the art. 

Since the impulse response of the channel is estimated for all the 
slots m of the received signal, it is possible to improve this estimate ft„ , by 
implementing methods known to those skilled in the art. 

The exemplary implementation of the method according to the 
invention, given by way of illustration, makes use of a Wiener estimation, 
known to those skilled in the art, the principles of which shall not be 
described in detail. The Wiener estimation method is described for example 
in A. Blanc-Lapierre and 8, Pichinbono, Fonctions Ateatoires (Random 
Functions), Masson 1981. Any other method of estimation that improves the 
channel estimate may be used without departing from the scope of the 
invention. 



The Wiener estimation of the channel 

In the Wiener estimation of the channel, the estimate h m is 
improved by means of a linear transform of H , . This transform is chosen so 
as to minimize the root-mean-square error between the true channel and the 
5 improved estimate. 

The improved estimate of the impulse response of the channel for 
a slot m takes account of the covariance matrix of the channel and the 
covariance matrix of the estimated channel and is written for example (9): 

10 TA- 1 ft. 

M 

S Where r-EFH H»l 

m r represents the covariance matrix of the channel : 1 Cl » m 

W with H m being the vector defined in the expression (8) which remains 

m 1 5 constant or practically constant throughout the duration of the slot, and 

" J A represents the covariance matrix of the estimated channel : A = E[H„H* ] 

JU The idea of the invention consists especially in estimating the 

S covariance of the estimation noise in order to deduct it from the covariance 

K matrix of the estimated channel A in order to deduce the covariance matrix of 

p 20 the channel r therefrom. Rom these two values, it is then possible to deduce 

™ the improved estimate of the Impulse response of the propagation channel. 

The noise estimation corresponds to the difference between the 

estimated channel H m and the real channel H m . This noise results 

especially from interference and other noise sources affecting the signal 

25 during the channel estimation procedure. 

To this end, the method uses the estimation of the statistics of the 

channel as described here below. 

Estimation of A, covariance matrix of the estimated channel 

This matrix can be estimated in a way known to those skilled in the art, 
30 assuming that, during a number M of slots, estimates of the channel 

have been accumulated for which the arrival time is have not changed or 
has hardly changed. By means of these estimates, with m varying from 1 to 
M and in using a principle known to those skilled in the art (the estimator 
chosen is the best unskewed estimator for example), it is possible to estimate 
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the matrix A containing the second-order statistics of the estimate of the 
channel. 

The matrix A is then estimated by (10) 

This estimator is consistent inasmuch as it converges on the true 
matrix Awhen the number of observations tends towards infinity. 

Estimation of the covarlance r of the channel 

tt is expressed in taking account of the estimated matrix (10) 




and 

of the estimated matrix of covariance of the noise (1 1) expressed in the form: 

N 0 P 

where R x is the empirical correlation matrix of the observations, 

Let L be the number of samples of the channel to be estimated 
and K the number of sensors, the matrix R x has a size KL*KL It may be 
considered to be a temporal matrix with a size L"L whose elements are 
space matrices sized K*K. 

On the ith line and the jth column, (i and j varying from 1 to L) of 
the temporal matrix, the space matrix (K*K) R X |, is estimated by the sum on n 
of the values x(n+i)x(n+j) H or, in an equivalent way, of the values x(n+i)x(n) H 
or again of the values x(n)x(n+j-i) H . 

The covariance matrix of the channel is estimated by (12) : 

r=— JM" ■ — r, 

R x being an estimator consisting of the matrix r, when the number of 
observations M tends towards infinity, this quantity converges on V. 

Estimation of the pulse response of the propagation channel 
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The puis© response of the channel H m is then estimated again 
from the estimated covariance matrix (12) and from the inverse of the 
covariance matrix of the estimated channel in applying, for example, the 
Wiener estimation method or any other method used to improve an estimate 
5 by means of its statistics. 

There is, for example, a known first estimated value H m for the 
propagation channel H OT , for example, the one used to estimate the 
covariance of the estimated channel and the Wiener estimator, namely the 
expression (9), is applied to deduce therefrom the new estimated value H m 
10 of the propagation channel H m (13). 

A summary of the different steps of the method according to the 
invention is given in figure 3, in the form of an algorithm comprising the 
following steps: 

f5 » Reception and sampling of the multisensor signal, (a), 

■ storage of the successive observation vectors, for example in 
the memory, (b) 

• estimation of the matrix R, , ( c) 

• estimation of the channel by correlation of the signal with the 
20 learning sequence , (d) 

• storage of the successive estimates u m of the channel h , 
(e) 

• estimation of the matrix A, (f) 

• estimation of the matrix r, from the estimates ( c) and (f), and 
25 by using the expression of the estimates of the noise in the 

form _I_r (g) 

■ improvement of the estimates of the channel by applying for 
example the Wiener method to the estimates (e), (f) and (g). 

The estimates are computed for example by means of a 
30 microprocessor programmed accordingly. This microprocessor may also be 
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adapted for the execution of the steps of the method according to the 
invention. 

The microprocessor may be connected with a processing device 
using the estimated values of the channel to implement a demodulation or a 
direction-finding operation as described very briefly here below. 

The method described here above in the context of downlinks can 
also be applied for uplinks, namely links from a mobile station to a base 
station, without departing from the framework of the invention. 

Figure 4 is a drawing of the demodulator in the case of an uplink. 

The mobile unit sends out a single control channel, containing a 
learning sequence for the channel estimation and signalling information, as 
well as one or more data channers. The data channels are all sent out at the 
same power, either on the in-phase channel or on the quadrature channel. 
The symbols sent on each channel are B-PSK (+1) (Binary Phase Shift 
Keying) symbols, 

The model chosen for the downlink can be adapted to the uplink in 
considering the control channel as the user being attended to, and the data 
channels as the other users of the cell. The other mobile units transmitting to 
the same base station are processed similarly to the other cells for the 
downlink. They are included in the total noise. 

The steps of the method described in the context of the downlink 
can be applied with these assumptions. 

Figure 5 shows an exemplary implementation of steps to reduce 
the rank of the matrices used during the method. 
For example, the steps run as follows: 

• reception and sampling of the signal, (a) 

+ estimation of the propagation channel (d,e) 
- estimation of the parameters of the channel (arrival time, power 
of the paths, etc) 

♦ building of the reduction matrix, 
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• reduction of the observations and of the estimates of the 
channel, in carrying out the steps described here above ( c) (f) 
(g) in the normal case, in using the reduction matrix , 

• Wiener estimation of the reduced channel. 

It is possible to improve both the efficiency and the complexity of 
this method by resorting to the reduction of the channel space. This 
technique consists in parametrizing the channel space by means of an 
orthonormal base constituted by a small number of vectors. 

Indeed, generally, in the context of radio channel transmission, 
owing to the specular nature of the channel (with a restricted number of paths 
in the channel), the matrix T is said to be a "deficient ranking" matrix. This 
means that the vectors H m are deduced from a linear combination of a small 
number of vectors u v u p forming a base of the channel space. If these 
vectors were to be known, it would then be enough to estimate the linear 

combination g m<1 g m ,p used to rebuild H m by : 

p 

p-l 

Let us group g ml g m , P in a vector g m and u p in a matrix 

U, We then have: 

H„ =U Sw 

The vector g m has a size (far) smaller than H m . It is therefore 
simpler to implement the Wiener estimator of g m than H m . It is possible 
indeed to estimate g,ti by: 

Hence, it is possible to estimate the co variance matrices A and S 
of and of 0 H Xm, where X m represents the vector of the observations 

[x m (0) H x m (L)T received at the slot m. It is then possible to estimate 
matrix £2 which is the covariance matrix of g m . 



and 
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For which a new estimate of the channel is deduced: 

In practice, the vectors m, .... u p are not known, but they can be 
6 estimated in several ways: 

- Eigen vectors associated with the greatest eigenvalues of the 
estimated matrix r. This technique must be implemented with a 
method for determining the number of eigen vectors to be 
chosen. 

io . Canonical vectors associated with the greatest values of the 

diagonal of T. This technique amounts to keeping only the 
coefficients of the channel whose power exceeds a threshold 
that is set beforehand. It is not necessary to estimate r 
entirely but only its diagonal, which is far simpler, 
is • Vectors associated with the shaping filter shifted by the delays 

of the channel. This technique can be implemented only if the 
delays of the different paths of the channel have first of all 
been determined. A standardization step must be planned. 
The rank reduction step can be implemented before the different 
steps of the method according to the invention, or again after the steps 
depending on the cases of application. 

The method according to the invention is used for example to 
demodulate signals or again to carry out direction-finding operations. 

Demodulation 

Figure 6 gives a view, in a graph of the binary error rate as a 
function of the signal-to-noise ratio, of the performance characteristics of the 
different estimates. 

The curve (I) corresponds to a classic estimator, obtained by the 
correlation method, the curve (II) being a curb obtained by implementing the 
steps of the method according to the invention and the curve (III) 
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corresponding to the results obtained by the additional use of a step for 
reducing the rank of the matrix. 

Many techniques of demodulation rely on knowledge of the 
propagation channel. The RAKE receiver, for example, estimates the 
transmitted symbol from the phase of the quantity: 
4(')=T^i>(Wo +nYc Q (nMk) ff x(lN 0 +n + k) 

When the propagation channel is not known, which is always the 
case in practice, an estimation of this channel is done. In figure 4, the 
performance characteristics of our estimator (and of its reduced rank version) 
are compared with that of the standard estimator for the RAKE receiver. 

The new estimator and, to an even greater extent, the reduced 
rank estimator, give major gains in terms of performance characteristics. 

Co-operative direction finding 

The invention can be applied especially to the improvement of 
direction finding measurements. 

Indeed, in direction-finding, the angles of arrival of the different 
propagation paths are determined by estimating the vector h at the different 
instants of arrival. The co-operative estimation of the arrival times and of the 
associated channels is done by means of a learning sequence. It is therefore 
possible to improve it by implementing the method according to the invention. 
The precision of estimation of the angles of arrival Is directly related to the 
quality of estimation of h. When the estimated channel is noise-infested, the 
measured angle is in a cone centered on the real angle. Improving the 
estimation of h reduces the difference between the measured angle and the 
real angle. This improves the resolution of the processing. 



